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Abstract 

The scaling limit of the two-dimensional self-avoiding walk (SAW) is believed to be 
given by the Schramm-Loewner evolution (SLE) with the parameter k equal to 8/3. The 
scaling limit of the SAW has a natural parameterization and SLE has a standard param- 
eterization using the half-plane capacity. These two parameterizations do not correspond 
with one another. To make the scaling limit of the SAW and SLE agree as parameterized 
curves, we must reparameterize one of them. We present Monte Carlo results that show 
that if we reparameterize the SAW using the half-plane capacity, then it agrees well with 
SLE with its standard parameterization. We then consider how to reparameterize SLE 
to make it agree with the SAW with its natural parameterization. We argue using Monte 
Carlo results that the so-called variation of the SLE curve with p = l/v = 4/3 provides 
a parameterization that corresponds to the natural parameterization of the SAW. 
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1 Introduction 



The scaling limit of the two dimensional self-avoiding walk (SAW) is believed to be the Schramm- 
Loewner evolution (SLE) with k = 8/3 ^3]- Previous Monte Carlo studies have compared the 
scaling limit of the SAW and SLEg/s as unparameterized curves [Hllini- These studies consid- 
ered random variables that did not depend on how the curve were parameterized and found 
excellent agreement between the distributions for the SAW and SLEg/s. The goal of this paper 
is to use Monte Carlo simulations to understand how the SAW and SLEg/a curves should be 
parameterized so that they agree as parameterized curves. 

We give a brief definition of the SAW. For a detailed treatment we refer the reader to |17j . 
The relationship of the SAW to SLEg/3 is the subject of fT^^. We restrict our attention to the 
SAW on a two-dimensional lattice restricted to the upper half plane. Fix a positive integer N. 
We consider all nearest neighbor walks which begin at the origin, have exactly N steps, remain 
in the upper half plane (except for the starting point) and never visit the same site more than 
once. There are a finite number of such walks, so we can put the uniform probability measure 
on this set of walks. We then attempt to construct the scaling limit by the following double 
limit. We first let — 00. For the case of the half-plane this limit was proved to exist in jl3j . 
but there is no general proof of its existence. We then take the lattice spacing to zero. This 
limit has not been proved to exist, but it is believed that it does and gives a measure on curves 
in the upper half plane which start at the origin and go to 00. The SAW curves are simple 
(do not intersect themselves) by definition. This does not insure that in the scaling limit the 
measure is supported on simple curves, but it is believed that it is. 

SLE was introduced in [22j and studied further in [21 . Expositions may be found in J2] 
and j21I. In the case of k = 8/3, SLE satisfies the restriction property and so belongs to the 
one parameter family of restriction measures studied by Lawler, Schramm and Werner |14j . 
They showed that the probabilities of certain events for SLEg/3 may be computed from certain 
conformal maps. In some cases they may be computed explicitly. Past Monte Carlo simulations 
have tested the equivalence of the scaling limit of the SAW and SLEg/3 by computing these 
probabilities for the SAW by simulation and comparing with the exact result for SLEg/3. Here 
are two examples. Let P be a point on the horizontal axis and consider the minimum distance of 
the SLEg/3 curve to the point. This is a random variable whose distribution may be computed 
explicitly for SLEg/3. For the second example consider a vertical line given hy x = c for some 
c 7^ 0. The SLEg/3 curve intersects it. Let Y be the distance from the horizontal axis to the 
lowest intersection. The distribution of Y may be explicitly computed for SLEg/3. Note that 
for both of these random variables the time parameterization of the curve plays no role. For 
these two examples and several others, excellent agreement was found between the SLEg/3 exact 
results and the simulation results for the scaling limit of the SAW P| IIU]. 

The SLE curve is usually parameterized so that its half-plane capacity grows linearly with 
the time parameter t. Before we take the scaling limit, the SAW has a natural parameterization 
given by the number of steps, or equivalently the distance along the curve since all steps are 
nearest neighbor. With a suitable scaling this should give a natural parameterization of the 



2 



scaling limit of the SAW. Equipped with these natural parameterizations, the scaling limit of 
the SAW and the SLEg/s parameterized curves are not the same. We illustrate this with a 
simulation in section [21 How should these curves be reparameterized so that they agree? 

There are two ways to ask this question. We can ask how to reparameterize the scaling limit 
of the SAW to make it agree with SLEg/s with its parameterization using half-plane capacity. 
The answer is obvious - one should use the half-plane capacity of the scaling limit of the SAW 
to parameterize it. We study this in sectional 

A less obvious question is to ask how we should reparameterize SLEg/a to make it agree with 
the scaling limit of the SAW with its natural parameterization. In section |3] we give simulation 
results that indicate that for the SAW with its natural parameterization, the p-variation of 
the walk exists and is a deterministic, linear function of time if we take p = Xjv = 4/3. In 
other words, one can recover the natural parameterization of the SAW by computing its 1/z/ 
variation. Some care is needed here since different definitions of the variation give different 
values. However, it appears they differ only by a multiplicative constant. In section El we use 
the 1/z/ variation of the SLE to reparameterize the SLE curve and compare this with the SAW 
with its natural parameterization. Section IHl discusses various computational aspects of our 
simulations of the SAW and the SLE. 



2 Need for a random reparameterization 

We review some facts about SLE. We refer the reader to jil2jj for more detail. We restrict our 
attention to chordal SLE in the upper half plane. It is defined by the differential equation 
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and the initial condition g^i^z^ = z. The dot denotes differentiation with respect to time t, Bt 
is a standard real- valued Brownian motion, and z belongs to the upper half plane H. For some 
initial points 2;, the solution only exists for a finite time interval. The set of initial points for 
which the solution no longer exists at time t is a random subset Kt of H. For k < 4, is 
just a simple curve 7(t), i.e., a curve that does not intersect itself. We use 7[0, t] to denote the 
image of the curve for the time interval [0,t]. So 7[0,t] = Ki. 

The parameterization of the curve 7 that results from eq. (0) corresponds to the half-plane 
capacity of the curve which is defined as follows. Since 7 is simple, H[\7[0, t] is simply connected. 
So there is a conformal map of this domain onto H. If we require the map to send 00 to itself, 
then this map has an expansion about 00 of the form 

00 

g{z) = cz ^ anZ'"^ (2) 

n=l 

The map is uniquely determined if we require c = 1 and = 0- The half plane capacity of 
7[0,t] is then the coefficient a^. We denote it by /icap(7[0, if:]). 
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The above definition applies to any simple curve 7. For the SLE curve we have 

hcap{-f[0,t]) = 2t (3) 

If we were to parameterize SLE so that 

hcap{-f[0,t\) =b{t) (4) 

for some differentiable, increasing function b{t), then the 2 in the Schramm-Loewner equation 
would be replaced by b(t). 

The scaling property of Brownian motion leads to a scaling property for SLE: 7[0, t] has the 
same distribution as -y/t 7[0, 1]. (Note that if we rescale a set by r, its half-plane capacity scales 
by a factor of r^.) Thus 

mm = ct (5) 

where c = _E[7(1)^]. 

The SAW has a natural parameterization. Let W{n) be the infinite SAW in the upper half 
plane on a lattice with unit lattice spacing. We extend W{t) to all t > by linearly interpolating 
between the integer times. The average distance from the origin to W{n) is believed to grow 
like n'^ with u = 3/4. The limit of taking the lattice spacing to zero can then be obtained by 
defining 

uj{t) = lim n-''W{nt) (6) 

n^oo 

This limit is believed to exist and give a probability measure on simple curves in the upper half 
plane starting at the origin, but this has not been proved. This definition of the scaling limit 
is analogous to Brownian motion as the scaling limit of an ordinary random walk. If we take 
W{n) to be an ordinary random walk in the plane and u = 1/2, then © would converge to 
two-dimensional Brownian motion. We refer to this parameterization of the scaling limit of the 
SAW as its natural parameterization. If this limit exists, then u}{t) has the same distribution 
as f^u^l). Thus 

E[u{tY]=bt^'' (7) 

where b = E[uj{l)'^]. 

The behavior of E['y{t)'^] and E[Lj{ty] show that to make the scaling limit of the SAW and 
SLEg/s agree we must at least do a deterministic reparameterization. If there is a reparame- 
terization, ^(t) = ■j{(j)(t)), of SLEg/s with a deterministic (pit) that makes 'j(t) and uit) agree 
in distribution, then we must have ii^[7(t)^] = E[uj{ty] and so the reparameterization must be 

m = iiat'n (8) 

with a = b/c. 

It is easy to see from simulations that this deterministic reparameterization does not work, 
i.e., 7(t) and uj{t) do not have the same distribution as parameterized curves. To show 
this we consider 7(1) and uj{l). If they have the same distribution, then 7(1)/£'[|7(1)|] and 
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Figure 1: The distribution, P{R < t), for the distance R from the origin to 7(1) (labeled by 
SLE) and to uj{l) (labeled by SAW). Both random variables have been rescaled so that they 
have mean one. 

ci;(l)/-E'[|u;(l)|] have the same distribution. For each of these two random points, we compute 
two random variables - the distance R from the origin to the random point and the polar angle 
of the random point. Our rescaling means that E[R] = 1 for both the SLE and the SAW. 

The distributions of R for the SLE and the SAW are shown in figure The distributions 
of are shown in figure El We do not show any error bars in the figures. The width of the 
error bars would be hard to see on this scale. Clearly the distributions are different. 

Throughout this paper our plots of the distributions of random variables are plots of their 
cumulative distribution functions rather than their densities. The cumulative distribution is 
the function that is actually computed in the simulation. Computing the density would require 
taking a numerical derivative. 

3 Reparameterizing the SAW 

In this section we reparameterize the scaling limit of the SAW so that it agrees with SLE with 
its usual parameterization. For t > we define Ct to be the random time for the scaling limit 
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Figure 2: The distribution, P(0 < t), for the polar angle 9 of the points 7(1) (labeled by 
SLE) and (labeled by SAW). 

of the SAW when 

hcap{uj[0,Ct]) = 2t (9) 

We define Lj(t) = uj(Ct) so that hcap{^Lu[0,t]) = 2t. Then Lj(t) and 7(t) should have the same 
distribution as parameterized curves. We test this by comparing the distributions of R and 
for 7(1) and = uj{Ci). As before, R is the distance from the random point to the origin 
and is the polar angle of the random point. 

If we plot the distributions of R for 7(1) and cD(l), the two curves are virtually indistin- 
guishable. So in figure El we plot the difference of these two distributions. The important 
thing to note about this figure is the scale on the vertical axis. The difference between the two 
distributions is typically on the order of a tenth of a percent. We caution the reader that in 
figure H the random variable R was rescaled so that its mean was 1 . This is not done in figure 
121 The mean of R here is slightly greater than 2. In particular, the sharp spikes in figure (21 just 
left of i? = 2 correspond to where the distribution is increasing sharply. In figure Q this sharp 
increase occurs just left of i? = 1. 

The difference of the distributions of for 7(1) and ct'(l) is shown injH Again, the most 
important feature of this plot is the scale on the vertical axis. These two distributions also 
differ by on the order of a tenth of a percent. 
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The error bars in figures El and 0] and subsequent figures are only statistical errors, i.e., the 
error arising from not running the Monte Carlo simulations forever. The error bars shown are 
two standard deviations. The standard deviation is estimated using the techniques of batched 
means. The simulations of the SAW and the SLE both contain systematic errors. These are 
discussed in section El 
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Figure 3: The difference of the distributions of R for the SAW parameterized using 

capacity, and 7(1), the SLE with its usual parameterization. 

To further test the equivalence of 7(t) and cD(t), we consider the covariances of these two 
processes. We think of these processes as taking values in rather than the complex plane 
and let 7j(t) and Ui{t) denote their coordinates with i = 1,2. The covariance of SLE is given 
by 

Q,{s,t) = E [yi{s)^M - E [7.(s)] E [^M (10) 
By scaling, the dependence of the last term on s and t is trivial: 

E [7.(.)] E [7,(t)] = v^i? [7.(1)] E [7,(1)] (11) 

So we will only study E [7j(s)7j(t)]. By scaling, 

E [l^{s)lM = si? [7.(1)7, (t/s)] =tEYi,{s/t)^,{l)] (12) 
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Figure 4: The difference of the distributions of B for u)(l), the SAW parameterized using 
capacity, and 7(1), the SLE with its usual parameterization. 

So it suffices to study this quantity with s = 1 and < t < 1. We can rewrite this as 

E [7.(1)7, W] = E hithM + E [(7.(1) - 7.(t))7,(t)] (13) 

By scaling, the first term is just Cij t where Cjj = E [74(1)7^(1)]. So we simulate just the second 
term. So let 

p.,(t)=i^[(7.(l)-7.(t))7,(t)] (14) 
For the SAW parameterized by capacity, we define 

a,,{t) = E[{uj,{l)-(b,{t))^M (15) 

The functions Pij{t) and crij{t) are shown in figure El for i = j = 1 and i = j = 2. (It is not 
hard to see that by symmetry these functions are zero when i 7^ j.) The functions agree so well 
that we have plotted individual points for (Tjj(t) and an interpolating curve for Pij{t) so that 
the two functions may be distinguished. It is much harder to simulate the SAW at a constant 
capacity than to simulate the SLE. Consequently the statistical errors for the SAW are larger. 
So we have only shown error bars for the SAW points. The difference between the covariances 
of SLE and the SAW parameterized by capacity is small, but appears to be greater than the 
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statistical error shown by the error bars in the figure. We emphasize again that these error bars 
do not include the systematic errors. Computing the capacity of a SAW is time consuming. 
This limits the SAW simulation to walks of modest length - one of the sources of systematic 
error. 




Figure 5: The covariance Pij{t) for SLE and the covariance o"ij(t) for the SAW using capacity 
as its parameterization. The top curve is i = j = 2, the bottom curve is i = j = 1. The SLE 
covariance is plotted with a curve and no error bars, while individual points with error bars are 
plotted for the SAW. 



4 2^- variation of the SAW 

We now consider the question of how to reparameterize SLEg/a so that it agrees with the 
scaling limit of the SAW with its natural parameterization. If one plots points on a SAW that 
are equally spaced in time, then the points appear to be equally spaced along the SAW path. 
In other words, the natural parameterization of the SAW corresponds to using the length of 
the walk as the parameter. Of course, this is only a heuristic statement since the SAW is not 
expected to have finite variation. However, we will see that the so called "p-variation" provides 
a way to make sense of the length along the SAW. 
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If X{s) is a stochastic process, then the p-variation is defined as follows. Let = tg < < 
'1^2 ' ' ' = t he a. sequence of partitions of [0,t]. Let n„ denote the nth partition and let 

||n„|| be the width of the largest subinterval of n„. We assume that ||n„|| goes to zero as n 
goes to infinity. We define 

var,{{X{s))o<s<t,n^) = J2\X{t]) - X{t]_,W (16) 

i=i 

With some abuse of the notation we will write farp((X(s))o<s<t, n„) as just varp{X{0, t), Yin). 
We define the p-variation to be the limit as the partition gets finer: 

varp{X{0,t)) = lim mrp(X(0,t),n„) (17) 

n— >oo 

We have not specified the nature of the convergence in the above definition, nor have we given 
any conditions on the sequence of partitions. In our simulations of the SAW and SLE we will 
use a sequence of uniform partitions for the n„. Changing the parameterization of the process 
is equivalent to changing the partitions used. As we will see, changing the parameterization 
can change the value of the p-variation. 

It is important to note that we do not take a supremum over all partitions in the definition of 
the ]9-variation. If one takes a supremum one obtains a different quantity which, unfortunately, 
is also called the p- variation in the literature. (It is sometimes called the strong p- variation.) 
For Brownian motion and p = 2 the variation we have defined is the quadratic variation studied 
by Levy jin]. If the sequence of partitions satisfies ||n„|| = o(l/log(n)), then the quadratic 
variation of Brownian motion converges with probability one to t jl]. But if the condition is 
relaxed to ||n„|| = 0(l/log(n)) this is not true j^]. For other values of p, the p-variation has 
been shown to exist for certain Gaussian processes including fractional Brownian motion and 
the local times associated with a symmetric stable Levy process. See the recent references 
[HI HH Uni and references therein. 

Now consider the SAW. If the scaling limit exists then \uj{tj) — a;(tj_i)| is of size {tj — tj^i)". 
So it is natural to consider the p- variation with p = We conjecture that 

mri/^(w(0,t)) = ct (18) 

for some constant c. If the 1/v variation exists for the scaling limit of the SAW, then it is 
trivial to see that 

E[varii^{u{Q,t))]= ct (19) 

for some constant c. So the content of the conjecture is that the 1/v variation exists and is 
non-random. For the remainder of the paper, p will be 1/v = 4/3, and so we will drop the 
subscript l/u on var. 

We support the conjecture with a simulation. We only consider a uniform partition with in- 
tervals of width At, and we take t = 1. In this case we denote uar(a;(0, 1), 11) by var{uj{0, 1), At). 
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The distribution of this random variable is shown in figure IHl for several values of At. Two dif- 
ferent sets of curves are shown in this figure. The curves on the right side of the figure are for 
the variation defined above. Note that the scale on the horizontal axis begins at 0.92, not 0. 
The figure clearly indicates that the distribution is converging to a step function as it should if 
the 1/z/ variation is constant. 

To study the convergence quantitatively, in figure [7| we plot the variance of the random 
variable f ar(ti;(0, 1), At) as a function of At. Three curves are shown in this figure. The 
variance of var{uj{0, 1), At) is the line of points labeled "var." In this log-log plot the data is 
very well fit by a line with slope 1. This corresponds to the variance being proportional to At. 
Note that this is what we would find if the process had independent, stationary increments. 
We do not expect the increments of the SAW to be independent, but this result suggests that 
the lengths of these increments are weakly correlated. 

The fractional variation defined above involves the parameterization of the SAW. A different 
parameterization could give a different value to this variation. (We will see an example of this 
later.) Another definition of the variation that does not depend on the parameterization is the 
following. Let At > 0. We define times tj as follows, to = 0. Given tj, we let tj+i be the first 
time after tj such that |a;(tj_|_i) — uj{ti)\ = (At)'^. The variation over the time interval [0,t] is 
then defined to be nAt, where n is the largest index with t„ < t. We denote this variation by 
f ar„o(ci;(0, t). At). The subscript no stands for "no parameterization," to emphasize that this 
definition does not depend on how the curve is parameterized. 

The curves on the left of figure El show the variation f ar„o(i^(0, 1), At) for several val- 
ues of At. The figure shows this variation is also converging to a constant. It also shows 
that the constant varno{uj[0,t]) is slightly less than the constant f ar(a;[0, t]). The variance of 
f arjio(i^[0, t]. At) is the line of points labeled "varno" in figure [7| This data is also very well fit 
by a line with slope 1. 

We now return to our first definition of the p- variation, eq. and consider what can 

happen if we use a different parameterization. We are particularly interested in what happens 
when we use the half-plane capacity to reparameterize the SAW. We denote the resulting 
variation by f arcap(a;(0, t), n„). We can numerically compute the capacity along a SAW and 
so we can simulate this random variable. We have simulated it with a uniform partition. 
Computing the capacity is difficult and this limits the lengths of the SAW's we can study. 
This in turn restricts the simulations to larger values of At. The distributions of this random 
variable for several choices of At are shown in figure |H| The first thing to note about this figure 
is the scale on the horizontal axis. These distributions have considerably larger variance than 
the curves in figure IHl with the same values of At. Nonetheless, it appears this this random 
variable is also converging to a constant as At goes to zero. We also see that varcap{uj{0, 1)) is 
slightly less than varno{^^{0, 1)). The variance of varcap{uj{0, 1), At) as a function of At is the 
set of points labeled "varcap" in figure [7| This data indicates the variance is converging to zero 
as At goes to zero, but it is not well fit by any line, at least for the range of At shown. 

For each of the three variations we have studied, the plots of their distributions for three 
values of At very nearly have a common point of intersection. This provides a simple way to 
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estimate the values of these variations in the hmit At ^ 0. We estimate 



var{u{0, 1)) 
varno{uj{0, 1)) 
var^p{uj{0, 1)) 



0.987 ±0.001 
0.972 ±0.001 
= 0.92 ±0.01 



(20) 



The error bars we have given are just guesses. As we discuss in section IHl there are systematic 
errors in the simulations which make it difficult to give reliable error bars. 
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Figure 6: The set of curves on the right are the distributions of var{uj{0, 1), At); those on the 
left are the distributions of f ar„o(co'(0, 1), At). 



5 Reparameterizing the SLE 

If the scaling limit of the SAW is indeed SLE, then the results of the previous section suggest 
that each of the l/u variations we considered should exist for SLE and give a parameterization 
which corresponds to the natural parameterization of the SAW. Of course, there is no way to 
compute the variation var{u![0,t]) for SLE since it requires knowing the parameterization we 
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Figure 7: The variance of the three different definitions of the p-variation as functions of At. 



are trying to find. We must either use f ar„o(t^[0, t]) or varcap{uj[0,t]). We use the former since 
its variance for a nonzero At is smaller. Let Vt be the random time on the SLE where 

varno{l[0,Vt]) = ct (21) 

where c = f ar„o(i^[0, 1]). We define ^{t) = •yiVt) so that varno{^[0,t]) = t>ar„o(u;[0, t]) Then 
j{t) and uj{t) should have the same distribution as parameterized curves. We test this by 
comparing the distributions of u^l) and 7(1) = 7(Vi). 

We must compute the constant c = varno{^[0, 1]) by simulation. This requires computing 
the variation f ar„o(ci;(0, 1), At) for several At and then estimating the limit as At —* 0. How- 
ever, when we compute the variation f ar„o (7(0, VJ]) we must also use a nonzero At and then 
attempt to take a limit as At 0. These extrapolations have some error in them. We attempt 
to minimize this error by the following trick. For a given At we define Vt by 

mr„,(7[0, Vt], At) = var^oHO, 1], A^) t (22) 

and then let 7(t) = 7(V(). In other words, for the constant c we use the estimate of c that comes 
from the SAW simulation using the same value of At that we use to compute the p-variation 
for the SLE. (Note that Vt now has some small dependence on At.) 
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Figure 8: The distribution of f arcap(u;(0, 1), At), the random variable that converges to the 
p-variation computed using the parameterization of the SAW by capacity. 

As before we compute the distributions of the random variables R and G. In figure IHl we 
show the difference between the distributions of R for a'(l) and 7(1). Two plots are shown 
corresponding to two different values of At. The difference is larger than that seen in figure El 
However, the figure also shows that the difference depends strongly on At. In figure^lwe show 
the difference between the distributions of B for uj{l) and 7(1). For the angle we see almost no 
dependence on the choice of At. 

Next we test the equivalence of 7(t) and uj{t) by comparing their covariances. As we saw 
before, the non-trivial part of these covariances is given by 

p.,(t)=E [(7.(1) -7.W)7,W] (23) 

and 

a,,{t) = E[Ml)-uj,{t))ujM (24) 

These functions are shown in figure ^2 for ^ = J = 1 and i = j = 2. The SAW and SLE 
covariances agree reasonably well. As in figure there is significant dependence on the At used 
to compute the p-variation of the SLE. For some values of At the agreement between the two 
covariances is not as good as that shown in the figure. 
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Figure 9: The difference of tlie distributions of R for co'(l), the SAW with its natural parame- 
terization, and 7(1), the SLE parameterized by the variation var no- 
li is tempting to compare the plots in figures and El The first thing to keep in mind 
is that in figure ^ t = 1 corresponds to heap = 2 while in figure ^2 ^ = 1 corresponds to 
var no = 1- Simulations show that the average capacity for a SAW with var no = 1 is close to 
0.5. Changing the capacity by a factor of 4 is the same as changing spatial scales by a factor of 
2, and so corresponds to a change in the covariance by a factor of 4. Indeed the vertical scales 
in the two figures differ by a factor of approximately 4. However, we should emphasize that the 
two plots use completely different parameterizations of the processes and so comparing them 
is not particularly meaningful. These two plots should be different even if they are rescaled to 
have the same vertical scale. 

6 Simulating SAW and SLE 

One method for simulating SLE is to approximate the conformal map gt by the composition 
of a large number of conformal maps which are known explicitly. A particular implementation 
was studied by R. Bauer [T]. This method is also briefly discussed for the radial case in pUj . 
Here we give a brief overview. A more detailed review may be found in |T. 

We consider the time interval [0, 1] and divide it into N subintervals. (We could use a 
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Figure 10: The difference of the distributions of B for uj{l), the SAW with its natural parame- 
terization, and 7(1), the SLE parameterized by the variation varno- 



uniform partition of [0, 1], but one obtains a more uniform distribution of points along the SLE 
curve by using a non- uniform partition of the time interval. See j7j for details.) The conformal 
map Qt can be written as the composition of conformal maps, each of which corresponds 
to the solution of the SLE equation over one of the small time intervals. Each of these N 
conformal maps is approximated by a simple conformal map. We use the conformal map that 
takes the half plane minus a slit onto the half plane. The parameters for these conformal maps 
(e.g., the length and angle of the slit) are chosen so that the conformal maps have the correct 
capacity and so that the driving function corresponding to the composition of these conformal 
maps approximates the original driving function a/k-B*. For example, this can be done so that 
the driving function of the approximation agrees with y/UBt at the times which are endpoints 
of the N subintervals. 

To compute a point on the SLE trace, one must apply 0{N) of these conformal maps. So 
the time to compute a single point is 0{N). If one wants to compute points on the SLE trace 
this will take a time 0{N'^). The paper [7] shows how to implement this method of simulating 
the SLE so that the time required to compute a single point is approximately 0{N^'^) rather 
than 0{N). 

Certain features of the SLE are relatively easy to study by simulation. For example, to 
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Figure 11: The covariance crij{t) for the SAW and the covariance Pij{t) for SLE using the 
variation varno for its parameterization. The top curve is i = j = 2, the bottom curve is 
i = j = 1. The SLE covariance is plotted with a curve and no error bars, while individual 
points with error bars are plotted for the SAW. 

compute the distribution of the SLE trace at a fixed time, one need only compute a single 
point on the SLE trace. One can use a relatively large value of N and generate a large number 
of samples. By contrast, computing the 1/z/ variation of the curve involves a double limit. One 
needs to let oo and then let the time interval. At, used to compute the variation go to 

zero. In practice this means that one must take a small At and then take A^ large enough that 
1/A^ is small compared to At. Thus the simulations that require computing the variation are 
among the most difficult. 

We simulate the SAW using the pivot algorithm jT7j. The particular implementation of the 
pivot algorithm we use is found in |8j. This algorithm is fast in two dimensions, and one can 
study quantities like the location of the SAW at a fixed time for walks with a million steps. 
The difficult part of the SAW simulations is computing the capacity of a walk. This can only 
be done for much shorter walks. 

We compute the capacity of a SAW using the zipper algorithm fTH. 1^ . We give a brief 
explanation of how the algorithm works in our context and refer the reader to [201 for more 
detail. Consider a curve uj{t) in the upper half plane which starts at the origin and ends at 
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P, e.g., a self-avoiding walk. To compute the capacity we need to find the conformal map g{z) 
that takes the half plane EI minus the curve onto H. It should be normalized so that g{oo) = oo 
and g'{oo) = oo. This determines the map up to a translation. The capacity is the coefficient 
of 1/2; in the expansion of g{z) about 00. It does not depend on the choice of translation. 

Let = zq, zi, Z2, - ■ ■ ,Zn = P he points along the curve. Let Ci denote the portion of the 
curve from Zi_i to Zi. Let gi{z) be the conformal map that takes EI\ Ci onto H. In addition to 
the above normalizations, we require gi{zi) = 0. Then gi maps C2 U ■ ■ • U C„ to a curve that 
starts at = gi{zi) and ends at gi{zn) and passes through the points gi{zi), i = 2, ■ ■ ■ ,n — 1. 
Now we let g2{z) be the conformal map that takes H \ (71 (C2) onto H. Then g2 o gi{z) takes 
C3 U ■ ■ ■ U C„ to a curve that starts at = (72 o S'l (^2) , ends at g2 o gi (zn) and passes through the 
points g2ogi[zi), i = 3, ■ ■ ■ ,n — l. We continue to remove one Ci at a time. gi{z) is the conformal 
map that takes HI \ gi_i o ■ • ■ o (^(2 o gi{C.i) onto H. The composition g{z) = gn ° ' ' ' ° 92 ° 9i{z) 
is then the desired conformal map. The capacity of g is the sum of the capacities of the gi. 

Note that even if the original curve is piece- wise linear, the segments of the form gi{Ci), 
92 ° 9i{C2) ■ ■ ■ are not. The key idea of the zipper algorithm is to approximate gi by the 
conformal map hi that takes EI \ onto EI where Di is a curve that starts at and ends at 
9i-i o ■ ■ ■ (72 o gi{zi). So Di has the same endpoints as gi^i o ■ ■ ■ o (72 o gi{Ci). The curve Di 
is chosen so that the map hi is relatively easy to compute. One possibility is to take D^ to 
be a line segment from to gi-i o • • • o 5(2 o gi{zi). Then there is an explicit expression for 
the inverse of hi and hi itself may be found by a Newton's method. (Some care is needed in 
implementing Newton's method 1201 •) Another choice is to take Di to be the arc of the circle 
that is perpendicular to the real axis and connects and gi^i o • ■ ■ o 5(2 o gi{zi). This has the 
advantage that hi only involves linear fractional transformations and a square root. 

Because of the need to compute gi_i o ■ ■ ■ o g20 gi [zi) for every point Zj, the zipper algorithm 
requires a time O(n^). The idea in that speeds up the SLE simulation may be applied here 
to dramatically reduce the time required by this algorithm. 

The curve Di starts and ends at the same points as gi-i o ■ • ■ o (72 o gi{Ci), but does not 
necessarily approximate it well. For a fixed curve one improves the approximation to the 
capacity by using more points zi, ■ ■ ■ , 2„ along the curve. We do not do this. Thus the capacity 
we compute is not exactly the capacity of the SAW. However, it is the exact capacity (up to 
numerical errors) of some curve that passes through the same lattice points as the SAW. We 
expect that this curve will have the same scaling limit as the SAW. 

As in any Monte Carlo simulation there is error because we do not run the simulation 
forever. We refer to this as "statistical error." It is relatively easy to estimate (we use batched 
means). There are also systematic errors in the simulations which we now discuss. 

For the SAW, the scaling limit is given by a double limit. First we should let the length 
of the walk go to infinity and then take the lattice spacing to zero. In practice we fix a large 
N and simulate SAW's with steps but only use the first A^' steps where A^' is significantly 
smaller than A^. The simulation is done with a unit lattice, but then we rescale the walks by 
a factor of {N')~'^. So the portion of the rescaled SAW that we use has a size of order 1 and 
approximates the scaling limit of the SAW with the natural parameterization running from 
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t = to t = 1. For the SAW simulations we will give the number of iterations of the Monte 
Carlo Markov Chain. This number is typically very large, but the samples generated are highly 
correlated. 

For the SLE, denotes the number of conformal maps used in the approximation. For the 
SLE simulations we will give the number of SLE's generated. These samples are independent. 

We end this section with a discussion of the parameters used in the various simulations. In 
figures □ and H the SAW is simulated using A^ = 1, 000, 000 and A^' = 200, 000. The simulation 
ran for 1 billion iterations. The SLE simulation used A^ = 100, 000, and 1 million samples were 
generated. These large numbers are possible because all we need to compute for the SLE is the 
location of its position at t = 1. 

In figures El and |3 the SLE simulation used is the same as that in figures ^ and El The 
SAW simulation requires computing the capacity of the walk. This limits the simulation to 
significantly shorter walks. We take A^ = 100, 000 and rescale the walk by a factor of . For 
this rescaled walk we compute the capacity along the walk up until heap = 0.1. (The number of 
steps at which this occurs is random but well below A^. The mean number of steps at which the 
capacity is 0.1 is about a third of the total number of steps.) We then rescale the SAW again by 
a factor of v^20 so that the random time we are finding is where the capacity is 2. As we have 
noted the samples generated by the SAW simulation are highly correlated. If we are studying 
an observable which is trivial to compute, we might as well compute the value of the observable 
at every time step. The samples will be highly correlated, but there is very little cost in terms 
of computation time. For observables that are not trivial to compute, e.g., observables that 
involve computing the capacity, if we computed the observable at every time step we would 
spend most of the computation time on computing the observable. For such observables we only 
compute it every s time steps. In this SAW simulation we took s = 100, 000. The simulation 
was run for 8 billion iterations, so 80, 000 samples of the observables were computed. 

For the first covariance plot, figure the SLE simulation is done with A^ = 10, 000. We 
compute the SLE covariance at 50 equally spaced times. 200,000 samples of the SLE were 
generated. For the SAW simulation, we take A^ = 100, 000, rescale the walk by a factor of A^"*^, 
and then compute the capacity along the walk up until heap = 0.1. We then rescale the SAW 
again by a factor of \/20 so that the random time we are finding is where the capacity is 2. 
Computing the capacity of the SAW is slow, so we only compute it every 100, 000 iterations of 
the Markov chain. We ran the chain for 8 billion iterations, so we generated 80, 000 samples of 
the SAW covariance. 

There are two simulations of the SAW in figure (HI In both of them A^ = 1, 000, 000 and 
we use A^' = 500, 000 and A^' = 1, 000, 000. For most observables, taking A^' at or near A^ is 
a bad idea - it will change the distribution of the observable. However, for the two variations 
shown in the figure we find no difference between the variations computed using A^' = 500, 000 
and A^' = 1,000,000. The curves shown use the latter value. Both simulations were run for 1 
billion iterations. 

In the simulation of the SAW used for figured we take A^ = 100, 000 and A^' = 10, 000. The 
relatively small value of A^' is because of the difficulty in computing the capacity of a SAW. We 
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run the Markov chain for 1 bilhon iterations, but only compute var^ap every 100, 000 iterations 
for a total of only 10, 000 samples. 

In figures IHl and ^1 the SAW simulation used is the same as that in figures Q and |2l For 
the SLE simulation we took N = 500, 000. However, we only computed every fifth point along 
the SLE, so we compute a total of 100, 000 points on the SLE. We generated 146,000 samples. 
This is the longest simulation in this paper. It required approximately 150 cpu-days. The 
time needed for this simulation is significantly reduced by the following trick. We compute the 
variation vavno up to the time when it is 1. This time on the SLE when this is attained is 
random, but occurs well before the end of the SLE we are computing. So we actually only need 
to compute some initial fraction of the 100, 000 points on the SLE. 

Finally, for figure EH the SAW simulation uses = 1,000,000 and A^' = 400,000. The 
SAW simulation was run for 1 billion iterations. For the SLE simulation, A^ = 250, 000 but 
only every fifth point was computed. Again, as discussed in the preceding paragraph we need 
only compute an initial portion of the SLE. 127,000 samples of the SLE were generated. 

7 Conclusions 

The main conclusion of this paper is that the p- variation with p = 1/z/ of the SLE trace 
provides a parameterization that corresponds to the natural parameterization of the SAW. 
Simulations show that the p-variation of the SAW is non-random. A trivial scaling argument 
then implies it is proportional to the natural parameterization of the SAW. Thus the SLE with 
the p-variation as its parameterization and the SAW with its natural parameterization should 
agree as parameterized curves. Two tests were done to check this agreement. The distributions 
of the random curves at a fixed time were compared, and the covariances of the two processes 
were compared. Good agreement was found. 

A secondary conclusion of this paper concerns using the half-plane capacity to reparame- 
terize the SAW. With this reparameterization it should agree with the SLE as parameterized 
curves. The same two tests of their equivalence as parameterized curves showed good agree- 
ment. 

For random fractal curves like the SAW or the loop-erased random walk the exponent v 
should be equal to the reciprocal of the Hausdorff dimension of the curve. Beffara [21 E] has 
proved that if 7 is the SLE trace with parameter k, then the Hausdorff dimension of 7 is 1 + k/8 
a.s. This suggests that the p-variation of SLE should exist for p = = 1 -|- /t/8. A natural 
question is to consider this variation for other discrete models and see if it exists and is non- 
random as it appears to be for the SAW. Preliminary simulations of the loop-erased random 
walk indicate that its 1/v variation exists and is non random for v = 4/5. This was one of the 
models considered by Schramm in his original paper j22j|. Lawler, Schramm and Werner J3] 
have proved that its scaling limit converges to SLE with n = 2. 
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